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SUMMARY 

An advanced parameter identification process was used to identify 
six d6f;ree-of-freedom rigid body stability and control derivatives from 
flight test measurements for a large "crane" type helicopter. Data was 
processed for a range of airspeeds and for pulse, step and doublet control 
inputs. The sensitivity of the identification results to processing 
parameters was investigated. The acceptability of the identified deri- 
vatives was examined using data regeneration and simulation, comparison 
with analytical results and using engineering judgment. 

The identified major derivatives agreed reasonably well with analytical 
results and intuitive engineering analysis. The non-major derivatives 
identified were inconsistent. The doublet input data provided superior 
results to the step and pulse input data. 

INTRODUCTION 

There are several existing techniques for determining values of 
stability and control derivatives for a fixed wing vehicle. They include 
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wind tunnel testing, semiemplrlcal methods such as U.S.A.F. DATCOM and 
Royal Aeronautical Society Data Sheets, and flight-test Identification 
methods (see Ref. 1). Wind tunnel testing can provide good estimates for 
static derivatives (e.g. , functions of a and 6) but has limitation for 
providing dynamic derivatives (e.g., associated with velocities or rates). 
Semlempirical methods which utilize aircraft geometry and aerodynamic 
characteristics have been used to provide preliminary design derivative 
values but are limited in accuracy and may not provide acceptable estimates 
for a new configuration. Semlempirical methods have not been developed 
for helicopter configurations and the same fixed wing limitations for wind 
tunnel testing exists for helicopters in addition to requiring more 
sophisticated models and rotor dynamic scaling. The third method, while 
existing for many years, has not enjoyed large success in fixed wing 
application and has only recently been employed for helicopter application 
(see Ref. 2, 3, 4). The method of determining stability derivatives from 
flight-test measurements of vehicle and control states is commonly known 
as "stability derivative extraction" and is a form of parameter identification. 
The general procedure for extracting stability derivatives from flight 
data is to form a small perturbation mathematical model of the aircraft 
equations of motion, provide measurements of the vehicle states, and solve 
for the coefficients in the aerodynamic portion of the model. In vector 
notation the governing sets of equations are: 

• 

X = AX + Bu 

where X contains the vehicle state variables, u the control system state 
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variables, A contains the stability derivatives, and B contains the control 
derivatives. The procedure Is to solve for the components of A and B given 
measurements of X, X, and u. 


SYMBOLS 

Values are given in both SI and U.S. Customary Units. The measurements 
and calculations were made in U.S. Customary Units. 

A lateral cyclic control angle, deg. 

Is 

A/s indicated airspeed, m/sec (knots) 

B longitudinal cyclic control angle, deg. 

Is 

f^ digital filter cutoff frequency, herf: 

f^ digital filter termination frequency, hertz 

g acceleration due to gravity, m/sec^ (ft/sec^) 

G process noise gain matrix 

IAS indicated airspeed, m/sec (knots) 

vehicle rolling moment of inertia, kg-ra^ (slug-ft^) 

vehicle product of inertia, kg-m^ (slug-ft^) 

lyy vehicle pitching moment of inertia, kg-m^ (slug-ft^) 

I^^ vehicle yawing moment of inertia, kg-m^ (slug-ft^) 

L total rolling moment, m-N (ft-lbs) 

M total pitching moment, m-N (ft-lbs) 

n white gausslan measurement noise vector 
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total yawing moment, m-N (ft- lbs) 
vehicle roll rate, rad/sec 
state variable covariance matrix 

cross covariance matrix 

vehicle pitch rate, rad/sec 

vehicle yaw rate, rad/sec 

X axis velocity component, m/sec (ft/sec) 

Y axis velocity component, m/sec (ft/sec) 

Z axis velocity component, m/sec (ft/sec) 
zero mean white gaussian process noise 
total longitudinal force, N(lbs) 
stability derivative vector 

total lateral force, N (lbs) 
total vertical force, N (lbs) 
vehicle angle of attack, rad. 
vehicle roll angle, deg. 
vehicle sideslip angle, rad. 
vehicle pitch attitude, deg. 
vehicle yaw attitude, deg. 

tail rotor collective pitch at blade root, deg. 

main rotor collective pitch at blade root, deg. 

designates control input perturbation 
derivative uncertainty, derivative units 
derivative variance, derivative units squared 
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Subscripts 


0 Refers to Initial trim condition. 


Superscripts 

A dot refers to first derivative with respect to time (e.g.. 


9u . 
“ 3t^ 


A bar refers to a vector quantity. 


Derivatives 

Derivatives are normalized by mass or inertia. Derivative units are 
determined by force or moment units divided by appropri.te state variable 
unit and mass or inertia, e.g. 1/sec. 


ft - lbs 

rad/sec - slug - ft^ 


ft slu 
lb sec 


1 = 1 / 


sec 


Derivative Notation 

A shorthand derivative notation is used to simplify derivative expressions 
as illustrated by the following example: 


In addition, four derivatives contain inertial terms and are defined as 
follows : 


3Y 


3Y 


3Z 


9X 

q 9q P 9p r 9r q 3q 
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Symbols Used on Computer Plots 


THETA COLL 
THETA TAIL 
1 SIGMA 


main rotor collective pitch angle, deg. 
tail rotor collective pitch angle, deg. 
derivative uncertainty (1 o) , derivative units 


BACKGROUND 

As previously stated, extraction of stability derivatives from flight- 
test measurements for fixed wing vehicles is not a new idea. Documented 
results (see Ref. 1) are found as early as 1925. Two of the predominant 
early methods were analog matching and least-squares curve-fitting tech- 
niques. Both methods have yielded some success in identifying primary 
derivatives and associated characteristic modes for uncoupled longitudinal 
or lateral-directional models. Applications of these early methods pro- 
vided a basis for development of advanced extraction techniques and have 
resulted in the identification of several inherent problems associated with 
estimation of stability and control derivatives from applications of flight- 
test measurement techniques. Some of these problems are associated with 
the assumptions necessary to establish a mathematical model of the vehicle 
and others are associated with the limitations and complications of in- 
flight measurement of state variables. The major problems involving the 
mathematical model are: (1) exclusion of terms which have non-negligible 

values, (2) lack of provisions for nonlinear aerodynamics, (3) exclusion 
of coupling (longitudinal and lateral-directional) effects, (4) accounting 
for random disturbances (wind gusts, etc.) and (5) lack of modeling higher 
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degrees of freedom. These errors are generally categorized as sources 
of process noise. Another inherent problem of the mathematical model is 
that there are more derivatives to be identified than there are Independent 
equations such that the solution (estimated derivatives) is nonunique. 

In addition, since the equations of motion are vrrltten about an equilibrium 
condition, the failure to precisely trim the vehicle introduces an error 
source. 

The problems associated with state variable measurements are: (1) 

Inaccuracy of physical transducers used to evaluate the variables, (2) 
errors introduced by filtering (recording, reducing, and preprocessing 
data), (3) measurement bias errors due to drift or calibration errors, 

(4) misalignment and instrumentation uncertainty, and (5) noise due to 
aircraft vibration, physical transducers, or recording system. 

In addition to these major categories of error sources, the choice of 
pilot input has been shown to be an important consideration. Thus, all 
three areas must be considered before successful derivative extraction from 
flight data is realized. 

Three types of methods for parameter identification have been developed 
over the years for aircraft application. They are: (1) equation error 

methods such as classical least squares which minimize process noise effects, 
(2) output error methods such as Newton-Raphson and Kalman filter which 
minimize measurement noise effects, and (3) advanced statistical methods 
such as Maximum Likelihood and Extended Kalman filtering which account for 
both process and measurement noise (see Ref. 1). Since helicopter flight 
data contain substantial measurement and process noise, successful extraction 
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requires application of one of the advanced methods or at least a combina- 
tion of high-quality data filters combined with an eqr-.tion error method. 

Since the helicopter stability derivative identification problem om- 
bines all of the previously mentioned problems, existing fixed wing 
analyses are not generally applicable. A parameter identification technique 
specifically designed for helicopter application (Ref. 2) has been applied 
to flight test data from a large "crane" type, articulated rotor helicopter. 
The method chosen, as described in Figure 1, consist of essentially four 
processes. Flight measured state variables are digitally filtered using 
a Graham high-quality, zero phase shift filter (see Ref. 5) with a 
very high roll-off rate. The purpose of the digital filter is to retain 
only rigid body frequency data. The digital filtered data is then pro- 
cessed twice through a Kalman filter, first to estimate measurement biases 
in the measured data set, and then to account for random discrepancies of 
related measurements and generally improve the quality of the data. The 
Kaima:. ♦'■'Iter can also reproduce states for which data is not available, 
provided related data is available. The Kalman filtered data is then used 
in a least squares algorithm to provide initial derivative and variance 
estimates for initializing the Extended Kalman filter (Bayesian Maximum 
Likelihood Estimator, Ref. 2) algorithm for parameter identification. The 
least squares estimated variances are known to be conservative (Ref. 6) 
and are adjusted upward. The parameter identification algorithm used 
determines the most probable estimates of the state variable and stability 
derivatives given the measured data. The parameters to be identified 
(stability derivatives) along with the aircraft states are combined to 
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form an augmented state; the resulting state equation Is nonlinear and 
an Extended Kalman filter Is used to approximate an optimal nonlinear 
filter. The resulting system equations are: 

X,(t) = f(X,, X ) + G, w (t) 

1 1 p 1 1 

X = 0 
P 


where i represents the maneuver number, represents the stability deriva- 
tive vector, G is the process noise given matrix and w is zero mean white 
gaussi^n process noise. The measurement equations can be represented by: 

^ white ga’ ^sian measurement noise and Z the 

measured state variable. The estimates for the state variable and stability 
derivatives are obtained from the following filter equations: 


X = f(X, X ) + P,,R"‘(Z - X) 
P ** 


X = P ,R~^(Z - X) 
P P * 


where P and P are the state variable covariance matrix and the crca" 

11 p, 

covariance matrix fcr the state and parameter vector and R is the measure- 
ment noise intensity matrix (see Ref. 2). The algorithm is programmed so 
that several data sets (maneuvers) can be processed simultaneously to 
minimize the effect of process noise. 

The validity of the identified derivatives is then evaluated using 
several procedures. The identified derivatives are combined with the 
chosen math model and state variable time histories are generated using 
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the control inputs from the Identification flight data (regenerated data). 

The Root Mean Square (RMS) error between the filtered flight data and the 
regenerated data is then computed. Time histories are also simulated 
using the identified derivatives and flight data control Inputs not used 
in the identification process (simulated data) and RMS errors are calcu- 
lated between the flight data and simulated data. The Extended Kalman filter 
algorithm estimates the derivative uncertainty (1 o) for each derivative 
estimate. The relative magnitudes of the estimated uncertainty and the 
time history of the uncertainty can be used to judge the results of the 
derivative identification process. In addition, the derivatives and 
chosen model are used to compute the vehicle characteristic modes (e.g. 
eigenvalues). Comparison of the computed eigenvalues and the identified 
derivatives are made with available analytical results and engineering 
judgment is used to assess the acceptability of the id> 'tified derivatives. 

Since the algorithm used has the capability to simultaneously process 
several maneuvers, data including all four pilot inputs can be used and a 
coupled six degree-of-freedom model can be identified. The aerodynamic 
portion of the model contains thirty-six stability derivatives and twenty- 
four control derivatives. The stability derivatives are first order Taylor 
series expansions for aerodynamic forces and moments for u, v, w, p, q, 
and r. The control derivatives are first order terms of longitudinal 
cyclic, lateral cyclic, tail rotor collective pitch and main rotor col- 
lective pitch. The equations representing the chosen model are given in 
Figure 2. 
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FLIGHT DATA 

A flight test of the large "crane" type, articulated rotor helicopter 
shown in Figure 3, was performed to gather data for application of the 
stability derivative identification method previously described. The 
vehicle characteristics are given in Table 1. The vehicle was instrumented 
to provide data for airspeed, angles of at; ck and sideslip, vehicle 
attitude (Euler angles), rates, linear and angular accelerations and pilot 
control inputs. Data were sampled eighty-three times each second and 
recorded on magnetic tape. The vehicle was flown at a nominal gross weight 
of 128998. N (29,000 lbs) with center of gravity at 8.33m (328 inches). 

The test variables were airspeed and pilot control input shape. The test 
procedure required trimming the vehicle at a selected airspeed and tlien 
independently apply sequential control inputs of longitudinal stick, lateral 
stick, pedal movement and collective stick displacement while returning 
the vehicle to trim between suc' assive control inputs. The matrix of test 
runs is shown in Table II. 


RESULTS AND DISCUSSION 
Procedure Variations 

Several variations to the general identification procedure described 
previously were conducted in order to assess derivative sensitivity to the 
identlf ication process. These included; variation of the digital filter 
cutoff and termination frequencies, adjustment of the leisl. squares esti- 
mated derivative variances and modifications to the least sc- res initial 
derivative estimates. 
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Digital Filter 

The ptirpose of digitally filtering the measured flight data is to 
remove high frequency data content to enhance rigid body derivative identi- 
fication. A previous application of the method employed in this study 
to another flight vehicle (Ref. 2) was performed with an arbitrary digital 
filter cutoff and termination frequency. Analytical predictions (see Ref. 7) 
of the rigid body characteristic mode frequencies for this vehicle are all 
less than 1 >.z. Flight data for this study was digitally filtered with 
cutoff frequencies of 1 hz, 2 hz, 3 hz, and 4 hz. The average R!HS error 
between the regenerated and filtered flight data, the major derivative 
estimates and uncertainties were tabulated and compared in order to assess 
the effect of the filter frequency. The major derivative estimates and 
associated uncertainties are shawn in Tables IV and V, respectively. The 
average derivative uncertainty did not appreciably change with filter 
frequency, but some major derivative values change substantially. At a 
cutoff frequency of 4 hz the major derivative estimates showed the largest 
variation indicating, that more consistent results would be obtained by 
processing data at lower than 4 hz cutoff frequencies. The digital filter 
cutoff and termination frequencies chosen determine the amount of high 
frequency rotor effects present in the filtered state variables and re- 
duces the overall measurement noise content. Figure 4 shows a time 
history of vertical acceleration measurement tc demonstrate the effect of 
the filter frequency on the measured state variables. The results shown 
in Tables III, IV and V indicate that the presence of high frequency data 
concent did not measureably degrade the estimation process below 4 hz 
digital filter cutoff frequency. 
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Least Squares Variance 


As previously stated the Extended Kalman filter algorithm used In 
this study requires estimates of initial derivative values and their vari- 
ances. The least squares algorithm provides the initial derivative esti- 
mates and an estimate of the derivative variances which is knovm to be 
conservative. Therefore, the normal procedure is to multiply all the 
least squares estimated variances by a multiplication factor. Data were 
processed for five Least Square Variance Factors (LSQVF) ranging from 
25 to 130. The RMS Errors between regenerated data and filtered data, 
major derivative estimates and uncertainties are tabulated in Tables VI 
through VIII, respectively. The RMS errors and derivative uncertainties 
suggest a good value for LSQVF is 100. Several derivative estimates, 

including Z and , were very sensitive to the LSQVF. 

MR 

Initial Derivative Estimate 

The sensitivity of the derivative Identification process to the initial 
derivative estimates was also investigated. The least squares derivative 
estimates were increased by fifty per cent and processed through the Extended 
Kalman filtering algorithm with a data set previously processed with the 
initial least square derivative estimates. The LSQVF for both runs was 100. 
The major derivative estimates for both data runs are shown in Table IX. 

Most of the final derivative estimates were within twenty-five per cent 
of the original final estimates. Several of the derivative estimates 
including the important derivative M^, showed a marked Increase. However, 
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the absolute value of M for this case was small and would not Indicate 

w 

a major chanee In predicted vehicle stability. 


Pilot Input Shape 


Flight data was obtained using three pilot input shapes. The data 
were processed in order to determine if there was any advantage in using 
a particular pilot input for identification among the three most common 
stability flight test Inputs. This analysis was not an attempt to 
determine an optimum input (see Ref. 8 and 9), but was intended to pro- 
vide an indication of the best data to be chosen from previous flight 
testing results. The three most common inputs used for stability and con- 
trol testing are pulse, step and doublet inputs. A typical input of each 
type as used in this study is shown in Figure 5. Flight data for a nominal 
60 KIAS case for each of the input shapes were processed, and the derivative 
estimates and uncertainties are shown in Tables X and XI. The RMS errors 
for regenerated and simulated data are shown in Table XII. The doublet 
input produced the lowest composite RMS errors for regeneration and simu- 
lation of other maneuvers. The derivative uncertainties were not a strong 
function of input shape. Examination of the derivative values shows that 
the predicted uncertainty was not a good indicator of accuracy since the 
derivatives estimated from data using step inputs are obviously in error, 
yet the predicted uncertainties are not high. For the step input run 
was estimated at -444.6 and AT 15.18. Both of these estimates are 


very high in comparison to all other results. The doublet input provides 
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the best results of the three shapes processed. This result is consistent 
with the intuitive reasoning that the doublet can provide good excitation 
while keeping the vehicle response within the linear assumptions of a 
small perturbation model. 


Airspeed Variations 

Data were obtained for indicated airspeeds of 0 m/sec (hover) , 51 m/sec 
(30 knots), 101 m/sec (60 knots), and 152 m/sec (90 knots), ''ariations of 
the major derivative estimates with airspeed are shown in Table XIII. Also 
shown in Table XIII are derivatives predicted by a digital flight simula- 
tion computer program (see Ref. 7) for the subject vehicle at the same 
flight conditions. The derivative values agree fairly well with the 
exceptions of L^, M^, and N^, and the control derivatives. The variation 
of identified derivatives with increasing airspeed also corresponds fairly 
well with the analytical computer program results. 

Final Data Sets 

Data were processed for the 101 m/sec (60 knots) flight condition Tor 
pulse and doublet inputs utilizing the results of the process sensitivity 
analysis. The LSQVF was set at 100 and the data was digitally filtered 
at 3 hz and 4 hz cutoff and termination frequencies, respectively. The 
results of this processing were analyzed using the following methods: 
examination of derivative values using engineering judgment and analytical 
predictions, RMS errors between regenerated data and filtered flight data 
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used in the identification process, BMS errors between simulated data and 
filtered flight data not used in the identification process, examination 
of estimated derivative variances and conversion characteristics and com- 
. arison of eigenvalues with analytical predictions. Since the rotor 
degrees-of-freedom are not modeled the estimated derivatives are a com- 
bination of the rigid body derivatives and average rotor contribution. 
Depending on the location of the rotor characteristic moues these deriva- 
tives may not be equal to the analytical quasi-static derivatives (see 
Ref. 10). This may be one cause for differences between the analytical 
and identified derivatives. Although engineering judgment cannot be used 
to assess the validity of each of the sixty derivative estimates, the 
i^ign and relative magnitudes of most of the major derivatives can be 
subjectively analyzed. The major derivative estimates for a pulse and 
doublet control input are tabulated in Table XIV, From engineering experi- 
ence (Ref. 11 and 12) it is known that should be negative since it is 
a direct function of the negative of drag coefficient and should be small, 

Z should * i negative and of moderate value since it is a strong function 
w 

cf the negative of lift curve slope. should be negative if the vehicle 
has longitudinal static stability, but since the rotor is destabilizing 
may be positive at low speeds until the horizontal tail becomes 
effective. M should be small and positive for an articulated rotor heli- 
copter a^ Mg should be negative since all helicopters are known to have 
pit''i damping. should be positive indicating directional static 
'■tabllity and should be negative since yaw damping is present. 
shculd be negative and fairly large since the large rotor provides good 
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roll damping. For the sign convention used the major control derivatives 

M- and should be negative and and L. should be positive, 

^is ®MR ®TR ^s 

with L. substantially greater than since the same control moment 

is available for rolling and pitching while there is a large difference 
in inertias. In addition, since the equations of motion for this investi- 
gation include inertial effects in the derivatives should be approxi- 
mately equal to minus forward velocity and should be approximately 
equal to forward velocity. Examination of Table XIV shows that the Identi- 
fied derivatives comply with all the intuitive tests except M^, which 
should be positive for both runs. Analytically predicted derivative 
estimates are shown along with the complete set of identified deri atives 
in Tables XV and XVI. The major identified derivatives compare reasonably 
well with the analytical derivatives but in general predict lower damping 
derivatives. In addition, the analytical program predicts a positive 
as would be expected. The non-major derivative values are not intui- 
tively known and the large variation in their estimates between data runs 
indicates a poor capability to predict these values. Table XVII presents 
RMS errors between regenerated data and simulated data, as compared to 
filtered flight data. Figures 6 and 7 contain the regenerated time 
histories for pulse and doublet Inputs, respectively. The simulated time 
histories for pulse input data simulated by doublet input identified 
derivatives and vice versa are shown in Figure 8 and 9. The regenerated 
and simulated time histories all diverge as a function of time, this is 
due to the unstable eigenvalues of the identified vehicle (see Ref. 2) 
combined with the presence of process noise. The time history error 
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equation has the same eigenvalues as the Identified vehicle and thus even 
If the derivative estimates are correct, the time histories will diverge. 

For the time histories generated in this study the regeneration and 
simulations were all forced to reinitialize at the filtered data values 
after each four seconds of flight data. As expected, the correlation 
between regenerated data and flight data used in the identification process 
is better than the correlation between simulated data and flight data not 
used in the identification process. The pulse input data estimates appear 
to be superior when examining only the regenerated results; however, the 
more stringent requirement of data simulation verifies the previous 
result that the doublet input provides the best identif iability. The 
derivative uncertainties (la) of the major derivatives are shown in Table 
XVIIl for both data runs. Predicted uncertainties for both data runs are 
relatively small with the doublet input run having slightly lower uncertain- 
ties. The high uncertainty associated with for the pulse run is con- 
sistent with engineering experience that should be positive. For the 
doublet input run the uncertainty of is very small and yet the estimate 
for was also negative. A low predicted uncertainty is therefore not 
considered to be a sufficient condition for derivative estimate acceptance. 
However, the uncertainty can be useful in determining whether or not 
enough data has been processed for derivative convergence or if the con- 
trols have sufficiently excited the vehicle to allow identifiability of 

control derivatives. An examination of Figure 10 shows that the M estimate 

u 

would be in doubt due to its convergence characteristics as compared to 
the estimate shown in the same figure. 
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A comparison o£ characteristic equation solution (eigenvalues) for 
the identified derivatives and analytical results is shown in Figure 11. 
The analytical and identified derivative eigenvalues had similar results 
consisting of three complex pairs and two aperiodic values. The results 
of Figure 11 are what would be expected from a classical fixed wing 
vehicle analysis with three stable complex pairs and one stable aperiodic 
mo'^e. Helicopters normally exhibit an unstable phugold oscillation. The 
identified doublet input derivatives provide the expected helicopter 
results while the pulse input data results in two unstable oscillations 
and the analytical program predicts a slightly stable aperiodic mode. 


CONCLUSIONS 

1. The digital filter cutoff and termination frequencies do not have a 
major effect on the derivative identification process used for this study 
below frequences of 4 hz and 5 hz, respectively. 

2. The Least Squares Variance Factor (LSQVF) had an appreciable effect 
on the derivative identification results. 

3. The major derivatives identified agreed fairly well with analytical 
predictions and complied with engineering judgment with the exception of 

4. The predicted uncertainty was not a reliable method for evaluating 
derivative acceptance but an examination of the time history of estimated 
uncertainty proved to be useful in determining the adequacy of data con- 
tent for acceptable derivative identification. 
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5. The derivative process used in this investigation was fairly sensitive 
to initial derivative estimates. 

6. The derivative process used did not produce consistent values for non- 
major (off diagonal) derivatives. 

7. Flight data regeneration using identified derivatives should be supple- 
mented by flight data simulation when analyzing identified derivative 
acceptance. 

8. Doublet control Inputs provide a better data set for identification 
than pulse or step inputs. 
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Fusgelage: 


TABLE I.- DESCBIPTiON OF TEST VEHICLE 


Gross Weight N(lbs) 

128998. 

(29000.) 

I Kg-m^ (slug-ft^) 

XX 

39/65. 

(29329.) 

I Kg-m^ (slug-ft^) 

yy 

20364. 

(150198.) 

1 Kg-m^ (slug-ft^) 
zz 

177611. 

(130999.) 

I Kg-m^ (slug-ft^) 
xz 

-18168. 

(-13400.)* 


Length m(ft) 


21.41 


(70.24) 


Main Rotor: Number of Blades 

Diameter m, (ft) 

Airfoil 
Chord cm (in) 

Total Blade Area m^(ft^) 

Disc Area m^ (ft^) 

Solidity Ratio 
Blade Twist (deg) 

Shaft Tilt-Longitudinal (deg) 
Shaft Tilt-Lateral (deg) 


6 

22.02 (72.26) 

NACA 0010.91 

60.07 (23.65) 

32.72 (352.20) 

378.30 (4072.00) 

.1021 
-7.2000 
3 
3 



Tail Rotor: Number of Blades 

4 


Diameter m(ft) 

4.88 

(16.00) 

Airfoil 

NACA 0012 


Chord cm (in) 

38.99 

(15.35) 

Total Blade Area m^(ft^) 

3.08 

(33.12) 

Disc Area m^(ft^) 

18.68 

(201.10) 

Solidity Ratio 

.198 


Blade Twist (deg) 

-8.000 







TABLE 1.- DESCRIPTION OF TEST VEHICLE 


Horizontal : 
Stabilizer 


Vertical : 
Fin 


(continued) 


Arean>2(ft2) 3.72 (40.00) 

Aspect Ratio 1.805 

Taper Ratio 1 . 0 

Incidence Angle (deg) 0. 

Dihedral Angle (deg) 0. 

Area ^^(ftZ) 1.55 (16.70) 

Aspect Ratio 5.34 

Taper Ratio .733 

Incidence Angle (deg) 0. 

Dihedral Angle (deg) 0. 


*Assumed zero for identification process 



TABLE II.- DATA RUN SUMMARY 


Flight 

Data Run 

Airspeed 

(IAS) 

Control Input Shape 

38A 

1 

i 

0 m/sec 
(hover*) 

Doublet 

38A 

4 

152 m/sec 
(90 Knots) 

Doublet 

38A 

5 

1 101 oi/sec 
(60 Knots) 

Puljd 

38A 

6 

101 m/sec 
(60 Knots) 

Step 

38B 

2 

51 m/sec 
(30 Knots) 

Doublet 

38B 

4 

1 

lOl m/sec 
(60 Knots) j 

Doublet 


*Out of ground effect. 


All IG Level Flight 


f 


I 


TABLE III.- REGENESATKW RMS ERRORS FOR 

FOUR DIGITAL FILTER FREQUENCIES 


Digital Filter 
f /f 

c t 

1 

1/2 1 

1 

4 

‘ 

2/3 

3/4 

4/5 

Attitudes 

.825 * 

.768 

1.370 

.750 

(deg) 

' 

i 




1 

1 Velocities 

; 

.234 

.236 



.237 

.230 

1 m/s (ft/s) 

(.767) 1 

(.773) 

(.777) 

(.754) 

Angular F.ates 

1 

.752 

.765 

.840 

.794 

(deg/s) 

4 





1 





i 

Linear Acc. 

.273 

.276 

.292 

.281 

m/s^ (ft/s*) 

(.895) 

(.907) 

(.959) 

(.922) 

A: ,ular Acc. 

2.170 

2.330 

2.040 

2.680 


(deg/s^) 






TABLE IV.- DERIVATIVE VALUES FOR POUR 
DIGITAL FILTER FREQUENCIES 


LSQVF = 100 


J 

Digital 

Filter 

X/. 

2/3 

3/4 

4/5 

X, 

1 

.04990 

-.043400 

-.095200 

-.01630 

* 

L 

.00140 

-.000035 

-.000549 

.00198 

-.00405 

-.002570 

-.003040 

-.00387 

n"" 

.00560 

.006640 

.006350 

.00552 

z 

-.32700 

-.428000 

-.224000 

-.57900 

m” 

w 

.00134 

-.000500 

-.000430 

.00434 

L 

-.65600 

-.755000 

-.607000 

-.28100 


.15700 

.152000 

.200000 

.38700 

Z 

92.10000 

89.300000 

80.930000 

110.80000 

q 

-.27100 

-.301000 

-.210000 

-.68700 

! L 

-1.15000 

-2.040000 

-1.790000 

-1.91500 

< 

-.65800 

-1.020000 

-.913000 

-1.00000 

Mb 

-.04090 

-.039000 

-.030000 

-.03760 

IS 


X 




.32400 

.333000 

.317000 

.27600 

IS 





®TR i 

.07300 

.058000 

.059000 

.06830 , 

^®MR 1 

-2.14000 

-2.800000 

-3.200000 

-3.05000 




TABLE V.- MAJOR DERIVATIVE UNCERTAINTY FOR 
FOUR DIGITAL FILTER FREQUENCIES 




Uncertainty 3 



Digital 

Filter 

c t 
*'u 

■ 1/2 

42. Y' 
; 40.0 

2/3 

. 

• 87.0 

. 34.8 

3/4 

13.1 

129.0 

i 

4/5 

1 

i i 

i 159.0 

j 26.1 

L 

40.2 

47.8 

11.0 

' 24.9 i 

: n'' 

11.4 

10.9 

6.2 

7.4 ' 

V 





Z 

9.6 

9.5 

3.4 

7.7 

\i'"' 
' V 

3.3 

18.0 

32.4 

iO.7 

L 

15.2 

15.6 

7.4 

20.9 


24.6 

16.2 

6.2 

8.0 

_ P 




! 

-7 

5.5 

6.5 

2.7 

! 6.0 


22.5 

18.6 

26.0 

' 11.3 

T 

4^ 

18.5 

12.8 

5.7 

6.9 

^’r 

14.0 

11.4 

6.3 

6.5 


7.6 

7.6 

5.1 

6.9 

s 





"a 

5.0 

4.3 

2.0 

3.0 

is 





N'. 





TR 

6.8 

7-0 

3.5 

4.5 

2^ 

”:-fR 

10.4 

. 10.9 

3.8 

7.4 


4 





‘ i 





TABLE VI.- REGENERATION RMS ERRORS FOR FIVE VALUES OF LSQVT 




LSQVF 

RMS 

r- • ’ r 

4 

25 1 

50 

75 * 100 j 150 

i _ _ i 1 i 

Attitudes 

(deg) 

1.89 

.938 

2.162 

1.370 

7 

* .892 i 

, ! 

t 

1 1 

■ Velocities 
, (ft/s) 

1.76 

.791 

2.375 

1.368 

.777 

! • ! 

. 798 

r 

i 

! A. Yel. 

1 (deg/s) 

1.18 

.852 

.840 

.848 

Linear Acc. 
(ft/s ) 

1.55 

1 

.973 

1.697 

.559 

‘ ■ ? 

.956 

Ang Acc. 
(de«»/s ) 

2. 84 

2.710 

i 

( 

1.253 j 
1 

2.984 ' 

! 

2.040 

1 

■ ■ '1 

2.7 40 

AVG 

1.844 

1 

2.117 ' 

1.197 1 

1.247 ■ 





TABLE VII.- DERIVATIVE 



ki. I 


FOR FIVE VALUES OF LSQVF 


75 

1 

1 

100 ' 

150 

118100 

-.090000 

-.101300 

-.000225 

-. 00 C 293 

-.000130 

-.005000 

i -. 005 f '>'' 

1 -.005160 

.005320 

.005330 

1 .005390 


t 


-.462000 

-.718000 

-.777000 

.001160 

.001180 

.000996 

-.515000 

-.509000 

-.543000 

.240000 

.278000 

.265000 

96.900000 

97.490000 

99.770000 

-.168000 

-.144000 

-.131000 

- 1.290000 

- 1.350000 

- 1 . 370 C '0 

-.679000 

-.695000 

-.664000 

-.042100 

-.037200 

-.036700 

.294000 

. 296000 

.295000 

.070500 

.070200 

.070100 

-.818000 

- 3 . 13 C 000 

- 3.260000 


4 




1 


I 



( 


'RODUCffilLITY OF THE 
.3INAL PAGE IS POOR 





TABLE IX.- EFFECT OF INITIAL DERIVATIVE 

ESTIMATE ON IDENTIFIED DERIVATIVES 


r 


Original 
Identified Derv. 

Identified Derv. 
with 50% Increase 
in Initial Estimate 

% 

Difference 

X 

u 

: 

1 -.095200 

I 

-.11080 

1 +16.4 

M 

u 

-.000549 

-.00044 

-19.9 

L 

V 

-.003040 

-.00620 

104.0 

N 

V 

-.006350 

.00515 

-18.9 

Z 

w 

-.224000 

1 

-.87100 

+289.0 

M 

w 

-.000430 

.00150 

+449.0 

L 

P 

-.607000 

-.50400 

-17.0 

N 

P 

. 200000 

.30000 

+5 .0 

Z 

q 

SO. 9 30000 

105.60000 i 

+30.5 

M 

q 

-.210000 

-.16400 

-21.9 

L 

r 

-1.790000 

-1.38000 

-22.9 

N 

r 

-.913000 

i 

-.73800 

-19.2 


”b 

IS 

-.030000 i 

-.03700 


+23.3 

^A 

IS 

.317000 j 

.29300 

1 

-7.6 

®TR 

.059000 1 

i 

.06990 


+18.5 

®MR 

! 

-3.200000 i 

-3.72000 


+16.3 


f 




TABLE X.- IDENTIFIED DERIVATIVE VALUES FOR 
THREE CONTROL HIPUT SHAPES 


Digital Filter - 1/2 LSQVF =62.5 


Input 

Pulse 

Doublet 

Step 

X 

-.030200 

-.25560 

-.98500 

u 

.000027 

-.00235 

.00086 

L 

-.001980 

-.00454 

-.00139 

V 

i 

.007180 

.00347 

.00356 

! z ! 

-.374000 

-.41600 

2.13500 

m” 

-.001170 

-.00048 

.00331 

. w , 

; L 

-.769000 

-.71700 

-.60800 

nP 

.117000 

.08080 

-.03050 

i P ! 

Z 

87.860000 

108.10000 

-444.60000 


-.239000 

-.75800 

-.94200 

q 




I L 

-1.740000 

-.10400 

-.62600 

1 

-.861CC0 

-.79800 

-.33200 

1 r 

♦ 





-.038300 

-.04390 

-.01750 

IS 




^A 

IS 

.335000 

.07700 

i 

. 18900 

TR 

.057300 

1 .05870 

1 

; .03860 

HR 

-2.720000 

1 -2.75900 

■ 15.18000 



TABLE XI.- DERIVATIVE DNCERTAIMTIES FOR 
THREE CONTROL INPUT SHAPES 


Uncertainty % 



j I 

' L j 12.8 67.9 10.7 

I 11.4 6.5 9.5 

r \ . . 






TABLE XII.- BMS REGENERATION AND SIMULATKB) 
ERRORS FOR THREE CONTROL INPUT 
SHAPES 


Flight Identified 

Data y from > Derivatives 


Input Shape 

38A-5 

Pulse 

38B-4 

Doublet 

38A-6 

Step 


Att. 

6.20 




4.910 


38A-5 

Vel. 

6.41 


1.858 


13.570 


Pulse 

Rates 
L. Acc. 

4.50 

4.37 

R 


s 

6.840 

16.680 

S 


A. Acc. 

4.67 




3.750 



Avg. 

5.07 


1.628 


8.247 



Att. 

15.10 




12.860 


38B-4 

Vel. 

14.78 


6.319 


85.100 


Doublet 

Rates 

8.63 

s 

2.641 

R 

11.810 

s 

Li* Acc* I 

10.43 


3.483 

45.370 



A. Acc. 

9.22 


2.766 


8.400 



Avg. 

1 

11.34 


3.766 


29.390 



Att. i 

13.94 


8.740 


41.600 


38A-6 

Vel. ! 

14.50 


9.190 


136.500 


Step 

Rates ' 

Le Acc« 1 

9.04 

10.23 

\ 

s . 

4.170 

6.075 

S 

29.500 

114.000 

R 


A. Acc. 1 

8.84 

i 

3.890 

1 

16.900 



Avg. 

11.14 

! 

6.260 

I 

> 

62.000 



R - Regeneration 
S - Simulation 
Digital Filter = 1/2 
LSQVF - 62.5 


RtlPRODUCIBILrrY OP THE 
Ot^i^AL PAGE IS POOE 














TABLE XIII.- EFFECT OF AIRSPEED ON DERIVATIVE VALUES 




Airspeed 

(IAS) 


Derivative 

0 m/sec 

52 m/sec 

101 m/sec 

152 m/sec 


(Hover) 

(30 Knots) 

(60 Knots) 

(90 Knots) 

i X 

-.07880 

.242000 

-.030200 

.053600 

u 

-.01890 

-.026000 

-.046400 

-.069900 

M 

.00805 

.004600 

.000270 

.001400 

u 

.00179 

.001210 

.000943 

.001320 

L 

-.04610 

.000287 

.001980 

-.005800 

V 

1 < 

-.03490 

-.028700 

-.020900 

-.014400 

! N 

-.00807 

.00217 

-.007180 

.004020 

1 V 

1 

.04100 

.03980 

.033500 

i .026200 

I 

1 z 

.12600 

-.10400 

' -.374000 

i .734000 

1 w 

i 

-.35900 

-.56400 

. 

1 -.800000 

1 

1 

1 

i 

-.949700 

j 

i M 

-.00249 

.00295 

I 

-.001170 

-.004200 

! 

-.00154 

-.00150 

-.001080 

j -.000808 

L 

-.88000 

-.48100 

-.769000 

-.361000 

P 

-2.19500 

-2.22400 

-1.992000 

-1.629000 

N 

-.24600 

-.29100 


-.117000 

: .008600 

j P 

1.10400 

.90200 

.481000 

-.001960 

: z 

-4.08000 

29.30000 

87.860000 

127.300000 ' 

i " 

1.35300 

51.39500 

100.650000 

j 150.020000 ; 

i M 

-.15500 

-1.03000 

-.239000 

-.149000 ■ 

1 ^ 

-.42700 

-.41400 

-.471000 

-.518000 1 

1 

L 

.41500 

.32400 

-1.740000 

1 

. -.134000 ! 

r 

.20400 

.29400 

.377000 

.613000 i 


.07020 

-.35700 

-.861000 

1 

-.794000 

N 

r 

-.36500 

-.47500 

; -.664000 

: -1.036000 

1 


-.03560 

! -.02870 

* -.030300 

j -.040800 


-.11240 

! -.11070 

‘ -.110900 

; -.114600 1 

1 

L. 

.37500 

.18100 


. 335000 

. 329000 

A 

IS 

.38100 

.37700 

: .375000 

. 382000 


.03630 

.03350 

.057300 

.074700 

®TR 

.09910 

. 10090 

.122600 

.143000 


-2.51000 

-3.23000 

-2.720000 

-2.900000 

®MR 

-5.64600 

-5.49300 

-6.494000 

-7.615000 


1st Item Identified derivative, 2nd item analytical estimates 


TABLE XIV.- MAJ* ’ IDENTIFIED DERIVATIVES 
FOR PULSE AND DOUBLET INPUTS 



Pulse Input 

Doublet Input 


38A-5 

38B-4 

X 

-.090000 

-1.23000 

u 

-.000293 

-.00321 

L 

-.005620 

-.00750 

V 

-.005330 

.00304 

Z 

w 

-.718000 

-.57500 

.001180 

-.00099 

L 

-.509000 

-.33100 

N^ 

P 

.278000 

.11300 

Z 

97.490000 

101.10000 

< 

-.144000 

-.72400 

L 

-1.350000 

.01660 

r 

-.695000 

-.78800 

Mr 

-.037200 

-.04400 

IS 



X. 

.296000 

.23200 

N 

®TR 

.070200 

.05800 

®MR 

-3.130000 

-2.86000 


Digital Filte’' * 3/4 
LSQVF * 100 
A/S » 


? 11 m/sec (60 Knots) 





TABLE (V.- PULSE INPUT IDENTIFIED DERIVATIVES COMPARED 
TO ANALYTICAL ESTIMATES (C-81) 



X 

Y 

Z 

L 

M 

N 

u 

-.09000 

-.02980 

.14700 

-.00965 

-.000293 

-.00263 


-.04320 

.07420 

-.08940 

-.01830 

.000943 

.01810 

V 

.02590 

.09030 

-.12500 

-.00562 

-.001070 

.00533 


-.00742 

-.24000 

-.05500 

-.02090 

.000946 

.03350 

w 

. 24000 

. 29000 

-.71800 

-.00860 

.001180 

-.00534 


.03960 

-.07250 

-.74600 

-.00293 

-.001080 

.01430 

p 

-2.86000 

2.97000 

37.00000 

-.50900 

-.229000 

.27800 


-.78700 

-5.68500 

-2.66200 

-1.99200 

.172000 

.48100 

q 

8.24600 

10.10000 

97.49000 

2.80000 

”1 

-.144000 

1.05000 


3.62700 

-4.71300 

97.92100 

-1.62600 

-.471000 

1.34300 

r 

-25.10000 

-95.59000 

-3.99000 

■ -~j 

-1.35000 

-.138000 

-.69500 


-.34300 

-96.33400 

3.16200 

.37700 

.037100 

- . 66400 

B 

1.79000 

-.32300 

1.91500 

. 10300 

-.037200 

.02060 

IS 

. 79900 

.02400 

1.97900 

-.04560 

-.111000 

.01060 

A 

-.24600 

-.00018 

-1.78400 

. 29600 

.008410 

.00873 

1 s 

.02580 

.58300 

.09240 

.37500 

-.003820 

.00322 


-.02410 

-.42800 

-.09930 

.04690 

.002950 

.07020 

TR 

-.02170 

-.41300 

-.00923 

-.07380 

-.004020 

.12300 


1.01000 

-.43100 

-3.13000 

02740 

.009500 

.03010 


.081900 

-.20500 

-6.04700 

.02520 

.020500 

. 08900 


Analytical Estimates are Listed Below Identi.ied Derivatives, Run 38A-5 
LSQVF = 100, Digital Filter - 3/4. 


I 


00 

I 

u 

u 

> c/i 

M Cs4 

HH 

§1 

M M 

CS H 
(d M 
C Cd 

Q >J 
U < 
M U 

(b h-4 

I- :-•. 

H>- 
2 hJ 
Id < 

0 2 

M < 

H S 
» c- 

& I— > 
2 :2 
l-t 

, a 

H U3 
Id Cd 

a < 
a a 

II 

1 


Id 


t 

! 

2 

-.00686 

.01810 

.00304 

.03350 

-.00648 

.01430 

.11300 1 

.48100 ■ 

1.12900 
1.34300 . 

-.78800 
-.66400 i 

1 

.02380 ! 

.01060 

.01900 

.00322 

.05800 
. 12300 

.02470 

.08900 

S 

-.003210 

.000943 

.001.740 

.000946 

r 

[ -.000990 
-.001080 

-.030. 00 
.172000 

-.724000 

-.471000 

. 

-.077800 

.037100 

-.044000 

-.111000 

—4 

. 010400 
-.003820 

.003560 

-.004020 

o o 
o 

lO 

! O 
C C 

• • 

r. 

-.07980 

-.01830 

-.07500 

-.02090 

-.02030 

-.00293 

-.33100 

-1.99200 

5.45000 

-1.62600 

.01660 

.37700 

o o 

O v£ 
<• m 

\0 >7 

o o 

• 

i ' 

— 

.23200 

.37500 

-.02450 
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o c 

o CN 
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■ c c 

• • 

i 1 
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-.45100 

-.08940 

.05910 
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101.10000 
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c ; T ' c. c: 

. -MO 

• » • 

: 1 

I- 

O C 

o a 
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t 1 

1 


.51500 

-.07420 

-.06260 

-.24000 

.03840 

-.07250 

15.47000 

-4.25500 

25.40000 

-4.71300 

-79.53000 
-97. 19400 

o o 
c c : 

O' <• 

O' CM 

o 

• *i 

1 

1 

-1.39000 
.58300 j 

1 

! o c 
} c o 

it 

‘ 5 sr 

I * * 

I I 

-.00128 

-.20500 

X 

1. . 

-1.23000 

-.04320 

-.03970 

-.00742 

-.36800 

.03960 

-16.14000 

-.78700 

21.97000 

2.19700 

-13.30000 

-.34300 

.09860 
. 79900 

i 

3.48000 

.02580 

• 

.24100 

-.02170 

.01400 1 
.08190 j 
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TASLE XVll.- REGEHERATED AND SIMDLATED 

ERRORS FOR PULSE AND DOUBLET INPUTS 
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TABLE XVIII.- MAJOR DERIVATIVE UNCERTAINTIES 
FOR PULSE AND DOUBLET INPUTS 
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Figiirf 2.- Ktiuations of motion. 
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Figure 6 
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figure 6.- Concluded. 
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Doublet input regenerated time histories 
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Figure 8.- Pulse input simulated time histories. 



PITCH ATT. DEGREES ROLL ATT. DEGREES YAW ATT. DEGREES 


1 



filtered flight data 


. . . . simulated responses 



igure 


o 


Doublet input simulated time histories 




N 

s 

M 


12 ^10 



TIME SEC 


rigure 10.- M and 
u 


‘'w t-inie histories. 




iiODHCIBlLlTV 

lOCHAL PA^GE ^ 



I 



Figure ll.- Eigenvalues for identXfi^' 
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